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We describe an iterative approach to computing long-time 
semiclassical dynamics in the presence of chaos, which elimi- 
nates the need for summing over an exponentially large num- 
ber of classical paths, and has good convergence properties 
even beyond the Heisenberg time. Long-time semiclassical 
properties can be compared with those of the full quantum 
system. The method is used to demonstrate semiclassical 
dynamical localization in one-dimensional classically diffusive 
systems, showing that interference between classical paths is 
a sufficient mechanism for limiting long-time phase space ex- 
ploration. 

Dynamical localization, the suppression of quantum 
phase space exploration in a system whose classical ana- 
logue is diffusive, is a remarkable example of non-ergodic 
behavior in quantizations of classically ergodic motion. 
Since its discovery almost two decades ago, the phe- 
nomenon has been discussed and observed in various nu- 
merical studies 0^], and also in experimental settings 
|3j. Formal connections with Anderson localization in 
disordered systems have been made Q . 

As expected from classical-quantum correspondence, 
diffusive quantum behavior is observed for times beyond 
the Ehrenfest time in systems with a diffusive classical 
limit. Localization then sets in (for dimensions d < 2) at 
a time scale whose dependence on the diffusion constant 
D and Planck's constant h can be understood by analyz- 
ing when interference between classically distinct paths 
begins to be statistically important. It is thus natural 
to ask whether phase interference between long classical 
paths alone is sufficient to produce localization, in the 
absence of "hard quantum" effects like diffraction and 
tunneling. Addressing this question has historically been 
made difficult by the exponential proliferation of paths 
with time in a chaotic system. For d = 1, for example, 
the localization time scales as D/h 2 , so exact semiclassi- 
cal calculations all the way to the localization time scale 
are in practice impossible to carry out for small values of 
ft, where the semiclassical approximation itself is likely 
to be valid. 

An attempt along these lines in Ref. || proved some- 
what inconclusive, although some preliminary evidence of 
anomalous long-time behavior was found at 21-27 time 
steps. This in itself was an impressive calculational feat, 
made possible by a symbolic dynamics and the piecewise- 
quadratic nature of the potential. In Ref. HfjJ, statistical 



properties of long periodic orbits were used to give plau- 
sibility arguments for semiclassical localization without 
performing explicit periodic orbit sums. In Ref. ||, the 
relationship was examined between the exact quantum 
propagator in a classically diffusive system, and the semi- 
classical one-step propagator. It is important to note, 
however, that iteration of a one-step semiclassical propa- 
gator does not bear much resemblance to long-time semi- 
classical dynamics, except insofar as both are (at least at 
short times) related to the quantum dynamics || . 

In this work we adopt a different approach, based on 
the idea, recently developed more fully in Ref. ||, that 
even though semiclassical propagation is not strictly mul- 
tiplicative, long-time semiclassical dynamics can in fact 
be well approximated by iteration of intermediate-time 
propagators, with controllable errors and well-defined 
convergence properties. In Ref. Q this approach was 
used to compute semiclassical dynamics past the Heisen- 
berg time Th (where individual eigenstates and eigen- 
values can be resolved), for a system with Th = 256, 
without exponential expenditure of computational effort. 
Good convergence properties with decreasing ft allowed 
for direct comparison between quantum and semiclassical 
stationary properties, such as long-time transport, spec- 
tra, and eigenstates. In the present work we do not di- 
rectly use any of the results of Ref. |j , but the interested 
reader is directed there for a more complete discussion of 
the underlying ideas. 

We begin by defining A t (i,j) to be the semiclassical 
propagator matrix taking quantum state j to quantum 
state i in time t (computed using the Gutzwiller-van 
Vleck semiclassical expression). Unlike the correspond- 
ing quantum propagator Ut(i,j), A t is not unitary, nor 
is it multiplicative, e.g. A t ^ (A t / 2 ) 2 . We can, however, 
easily estimate the deviation from exact multiplicativity 
of the semiclassical propagator, at least in the caustic- 
free case. We define the natural basis-independent L 2 
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where the normalization ensures that the norm of a uni- 
tary operator is one. We can then write 



\\A t -(A t/2 )Y = 0(h a ) 



(2) 
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where for smooth dynamics, we have the exponent 
(^smooth — 2. This can be seen by noting that A t 
is exactly given by combining two A t / 2 propagators, 
as long as the intermediate integration at time t/2 



1 



is performed by stationary phase. The relative error 
between performing the intermediate integral exactly 
[{A/2) 2 (x,y) = J dzA t / 2 (x, z)A t / 2 (z,y)] and by station- 
ary phase [A t (x,y) = J ap dzA t / 2 (x,z)A t /2(z,y)] scales as 
h (this being the order of the subleading term in the sta- 
tionary phase expansion); thus a smoo th = 2 |1. 

In the case of a discontinuity in the underlying clas- 
sical dynamics, which is to be considered in the present 
work, the situation is rather different. However, the gen- 
eral approach applies also to this case. In expanding 
around each of the stationary paths that contribute to 
At(x,y), the region around the stationary phase point 

1/2 

where the phase is slowly varying scales as U 1 . Thus, 
for a stationary intermediate point y within 0(Ti 1 ^ 2 ) of a 
discontinuity in the potential (or in its first derivative), 
the exact integral gets cut off within the region of gaus- 
sian integration, and the relative error between the full 
integral and the stationary phase approximation is of or- 
der unity. It is easy to see that for small Ti this diffractive 
effect dominates the effect of the subleading terms in the 
expansion (which as we saw lead to a smoo th = 2), and 
results in 



^disconti 



1/2. 



(3) 



For long times t, many classical paths must be summed 
over to obtain the semiclassical propagator; here we are 
of course assuming that errors in the sum over paths 
add no more coherently than the actual contributions 
themselves. 

The next step is to extend Eq. || to the more general 
form 



At - (A t/M ) M f = 0(Mh a ) = O 




th a 



(4) 



which follows from assuming the successive errors in re- 
placing M — 1 stationary phase integrals by exact ones 
to add incoherently. The assumption of incoherent addi- 
tion of errors breaks down for very large M (specifically, 
for M greater than h~ 1 , the Heisenberg time measured 
in units of the shortest periodic orbit fe[). However, the 
higher-order corrections in M will not be relevant for our 
purposes. 

Eq. ^ allows successive controlled approximations to 
be computed to the exact semiclassical dynamics by tak- 
ing t/M ^> 1 (this produces values much closer to the 
semiclassical than to the quantum results). Of course, 
taking M — ► oo (t/M <C 1) in the expression (A 4 /m) M > 
we instead recover the quantum propagator, as in the 
Feynman path integral formalism. The intermediate 
case t/M ~ 1 (as in the Bogomolny surface of section 
approach [ [Tcfl ) produces a long-time dynamics which is 
strictly speaking neither quantum nor semiclassical, and 
provides an interpolation between the two worlds. 

The scaling properties of the iterative semiclassical ap- 
proximation with h, time t, and "quantization time" 



Tq = t/M, as expressed in Eq. ^ above, hold even for 
times t beyond the Heisenberg time of the system. These 
scaling properties, based on power-counting arguments, 
have been extensively tested numerically in |9| . One qual- 
ification is that while for a — 2 and Tjj ~ h 1 the ap- 
proximation at the Heisenberg time using fixed Tq gets 
better and better as h — * 0, in the present situation we 
have a = 1/2 and Tjj ~ h~ , so as Ti gets small we need 
larger Tq to preserve the accuracy of the approximation. 

We are now ready to apply the above outlined formal- 
ism to the case at hand: dynamical localization in one- 
dimensional systems. We consider a kicked map 0Jll|] on 
a cylindrical phase space < q < 2n, — oo < p < oo, 
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p-V'(q) 

q + p mod 2n , 



(5) 



with kick potential 
V(q) 



\K{q 



Bcosq (0 < q < 2vr) (6) 



turned on momentarily once every time step. Locally 
the dynamics looks everywhere like an inverted harmonic 
oscillator with a sinusoidal perturbation (as long as B < 
K), except for a discontinuity in the impulse at q = 0. 
The classical motion is completely chaotic, and diffusive 
in p, 



(iP -Po) 2 )classical = Dt , 

with diffusion constant 



D = {{V?) = ~ 
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The quantization of Eq. || is straightforward [p],pr|. 
Choosing periodic boundary conditions in g-space, we 
have a momentum basis given by p n = nK, n = 
— oo . . . oo. The dynamics is given by a unitary one-step 
propagator 



U 



e -tp 2 /2h e -lV(q)/h _ 



(9) 



Because the quantum dynamics (as well as the classi- 
cal) is symmetric under parity [p — * —p, q — * 2ir — q], 
we will in what follows focus only on the even sector 
b)even = (\p) + \ ~p))/V^, P > 0. This eliminates the 
problem of tunneling between positive and negative mo- 
menta. 

The semiclassical dynamics (in the absence of caus- 
tics and Maslov phases, which are conveniently avoided 
by taking B < K in Eq. ||) is given by the standard 
Gutzwiller-van Vleck propagator 



A sc (p',p,t) 
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det 
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x exp 



iSj{p,p',t) 
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where Sj is the action for classical path j taking p to 
p' , and the determinant is the corresponding classical 
probability density. We can now use Eq. [l(] to evalu- 
ate the semiclassical propagator matrices A t q for various 
"quantization times" Tq M, and then iterate to obtain 
A t ~ (At q )'/ Tq . As Tq — > t, we obtain the exact semi- 
classical behavior. In general, though, we only need to 
take Tq large enough to obtain the desired level of con- 
vergence to the true long-time behavior. The behavior 
of the iterative semiclassical limit as Tq becomes large 
can be compared with the quantum dynamics as given 
by U t = U* (Eq. §). 

The numerical results of this localization study are pre- 
sented in Figs. 0, ||. In Fig. [l| we choose a piecewise linear 
map, with parameters K = 1.073, B = 0.0 in the kick 
potential of Eq. [|. In Fig. |[ a sinusoidal term is added 
to the potential: K = 1.073, B = 0.52. In both cases, 
the quantum and semiclassical calculations are performed 
with Ti — 0.293 (note that this takes us well into the 
semiclassical regime: the relevant expansion parameter 
is Ti/(2w) 2 , since (2tt) 2 is the area of a unit cell in phase 
space). We can now compute ((n — no) 2 }, the spread 
in momentum space in units of h, as a function of time 

(n =p/h). 

More explicitly, given a propagator G t , whether quan- 
tum or semiclassical, we define 



((n-n ) 2 ) = 
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The average over j is performed over initial momenta far 
from and also far from the edge of the numerical lattice. 

The classical diffusion result given by Eqs. 0, || ap- 
pears on the log-log plot in Figs. |l|, || as a straight line 
of slope one. The full quantum calculation is seen in 
the dashed curve, which in each case is seen to turn 
over and approach a constant after the localization time 
Tioc = D/h 2 . This theoretically expected value of T\ oc 
(which is also equal to the predicted RMS spread in mo- 
mentum at infinite time in units of h, the square root of 
the quantity plotted in Figs. 0J|) is given by T\ oc = 44.1 
in Fig. and T ioc = 32.7 in Fig. [| All these predic- 
tions are in reasonable agreement with the full quantum 
numerics. 

We now proceed to the semiclassical analysis. The 
dotted curve in each of the two Figures represents the 
momentum spreading given by iteration of the one-step 
semiclassical propagator A\ . This does not closely follow 
the full quantum (or, as we shall soon see, the exact 
semiclassical result), and the behavior of this quantity in 
Fig. H is particularly erratic. 

Now, guided by Eq. [|, we look for convergence to 
the exact semiclassical answer as the quantization time 
Tq = t/M is taken to be much greater than one. Specif- 
ically, in Fig. |l| we plot as solid lines the calculations 
with Tq = 7, 8, 9 (i.e. we take successive approximations 
A t (A 7 )'/ 7 , etc.). We see that the agreement between 



the three calculations is very good, strongly suggesting 
that convergence has been achieved. The semiclassical 
calculation begins to deviate from the quantum sometime 
around Xi oc ; nevertheless it does very clearly localize at 
a well defined momentum spread somewhat larger than 
that given by the quantum calculation. (Some differ- 
ence in the details of the end of classical diffusion for the 
quantum and semiclassical calculations is not surprising. 
One should note here for example that the discontinuity 
in the kick potential will have a diffractive effect on the 
quantum dynamics, one which will not be present in the 
semiclassical approximation.) 

In Fig. ||, the Tq — > oo convergence to the exact long- 
time semiclassical dynamics is found to be somewhat 
slower (calculations with Tq = 8,9,10 are plotted as 
solid lines). Nevertheless, up until t ~ 200 w 6Xi oc , 
the three curves are in very good agreement, with un- 
certainty small not only compared to their common de- 
viation from classical diffusion (straight line), but also 
compared to their common distance from the quantum 
curve (dashed). The evidence for localization is very clear 
in this case also, as is the failure of the one-step itera- 
tive approximation (dotted curve) to reproduce long-time 
semiclassical behavior. 

As an additional test of semiclassical localization at 
very long times, we consider the eigenstates of the suc- 
cessive propagators At q as Tq — > oo. In the absence 
of interference effects (i.e. considering each At q simply 
as a band random matrix of band width ~JDTq/K), we 
would expect the momentum spread (Sp) 2 of the typical 
eigenstate to increase linearly with Tq. In fact, however, 
phase interference between classical paths turns out to 
be very important indeed, and the average RMS width 
of the eigenstates of At q is found not to increase signifi- 
cantly with Tq once Tq > 1. In Fig. |, the mean RMS 
width (in units of %) of At q eigenstates centered well 
away from the edge of the numerical lattice is plotted vs. 
the quantization time Tq . Parameters are the same as in 
the previous two figures: plusses are used for the B = 0.0 
case (corresponding to Fig. [|) and squares for B = 0.52 
(as in Fig. g). The result using the Tq = 1 one-step it- 
erated semiclassical propagator is also plotted, as is the 
quantum momentum spread at Tq — 0. We see localiza- 
tion at large Tq for both sets of parameters; in each case, 
the localization length is somewhat larger semiclassically 
as compared to the quantum calculation. The result ob- 
tained using Tq = 1 is intermediate between semiclassical 
and quantum in both cases. 

We can test the scaling of the semiclassical localiza- 
tion length with the diffusion constant by comparing the 
results for B = 0.0 and B = 0.52. The ratio of semi- 
classical localization lengths for these two parameters is 
found to be 1.32 (using Tq = 6); analytically we predict 
1.35. 

The behavior of the semiclassical localization with h 
has also been investigated, and we find that the local- 
ization length increases roughly in accordance with the 
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theoretical predictions of dynamical localization theory. 
However, the semiclassical localization length does grow 
somewhat more slowly with \/U than the quantum local- 
ization length, apparently leading to a convergence be- 
tween these quantities at small h. Thus, for B = 0.52 and 
using Tq — 6 we find a ratio of 1.66 between the long-time 
semiclassical and quantum localization lengths; reducing 
Ti by a factor of 2 causes this ratio to drop to 1.35. Unfor- 
tunately we were not able to investigate extremely small 
H due to computer limitations. 

Thus, using an iterative approach to long-time semi- 
classical calculations, we have been able to see explicitly 
semiclassical localization in classically diffusive systems 
at small Ti. We can now say definitively that although 
details of long-time quantum dynamics are affected by 
diffraction and tunneling corrections, the essence of the 
localization phenomenon is indeed contained in the in- 
terference among long classical paths. 
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FIG. 1. Momentum spread ((n — no) 2 } (as defined by 
Eq. [H]), as a function of time for a kicked system with kick 
potential parameters K — 1.073, B — 0.0 (Eq. ^|). Classical 
diffusion (Eq. [?]) appears as a straight line; the quantum calcu- 
lation, represented by a dashed curve, shows dynamical local- 
ization at time scale T\ oc = 44.1. Successive approximations 
to the long-time semiclassical propagator (using Tq = 7, 8, 9) 
are drawn as solid curves. The one-step iterated semiclassical 
propagator (Tq — 1) produces the dotted curve. 

FIG. 2. Same as previous Figure, with sinusoidal pertur- 
bation in the potential, B — 0.52. Here the expected value of 
Tioc is 32.7. Approximations to long-time semiclassical prop- 
agation using Tq = 8, 9, 10 are shown as solid curves. 

FIG. 3. Average RMS eigenstate width for succes- 
sive approximations to long-time semiclassical propagation, 
Tq = 3 ... 9. The average eigenstate width for the one-step 
iterated propagator (Tq = 1) is also displayed, as is the quan- 
tum result (Tq = 0). Plusses represent the case B — 0.0; 
squares represent B = 0.52. 
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